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ABSTRACT 



It has been suggested that the mass-to-light (M/L) ratio of gravitationally 
clustering objects is scale-independent on scales beyond galaxy clusters, and may also 
be independent of the mass of the objects. In this paper, we show that the scale 
behavior of M/L ratio is closely related to the scaling of cosmic structures larger than 
clusters. The scale dependence of the M/L ratio can be determined by comparing 
the observed scaling of richness function (RF) of multi-scale identified objects with 
the model-predicted scaling of mass function (MF) of large scale structures. Using 
the multi-scale identified clusters from IRAS 1.2 Jy galaxy survey, we have made 
comparisons of the observed RF scaling of IRAS r c ;-clusters with the MF scalings 
given by simulations of three popular models SCDM, LCDM and OCDM. We find 
that, the M/L ratio basically is scale-independent from the Abell radius up to about 
24 /i _1 Mpc, while it seems to show a slight, but systematical, increase over this scale 
range. This result is weakly dependent on the cosmological parameters. 



Subject headings: cosmology: theory - galaxies: fundamental parameters - large-scale 
structure of universe 



1. Introduction 

The scale-dependence of mass-to-light M/L ratio is always an important topic of cosmology. 
The M/L ratio on scales of galaxies, double galaxies, groups of galaxies, and galaxy clusters 
provided the first approach of evaluating the mean mass density of the universe (e.g. Faber & 
Gallagher 1979). A most striking result in this study is that M/L ratio increases with the scale 
of the system considered (e.g. Rubin 1993). This has inspired the searching for dark matters in 
objects on scales larger than clusters. Yet, the M/L ratio appears to be flat on the scales of rich 
clusters and larger. This result has been employed to propose that the M/L ratio on scales larger 
than clusters asymptotically reaches a constant, which is approximately equal to the median M/L 
ratio of rich clusters, i.e. M/L ~ 300 ± 100 h (M/L) (Bahcall, Lubin & Dorman 1995). This 
asymptotic M/L ratio yields O m ~ 0.2 ± 0.07, suggesting a low density universe (Bahcall & Fan 
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1998). That is, the behavior of M/L ratio on large scales is directly related to the cosmological 
parameter U m . 

The evidence for a flat M/L on large scale is still weak. The observed Virgo cluster infall 
motion seems to be the only evidence of the flatness on the scale of the local supercfuster (~ 10 
h^ 1 Mpc) (e.g. Yahil, Sandage & Tammann, 1980). However, a direct measure of the mass-to-light 
ratio of the Corona Borealis supercluster gives M/L ~ 560 h{M/L) & . Namely, the M/L ratio on 
scale ~ 20 h~ l Mpc is higher than that of rich clusters by a factor of 2 (Small et al. 1998). On 
the other hand, galaxy clusters with Abell radius generally are dominated by elliptical galaxies 
at their centers, while larger structures have a higher fraction of spiral galaxies. Therefore, it 
is uncertain if the morphological segregation leads to a variation of the M/L ratio of structures 
larger than clusters. 

In this Letter, we try to approach this problem by analyzing the observed scaling of clustering 
on scales beyond clusters. We have found recently that the richness function (RF) and the mass 
function (MF) of objects identified by multi-resolution analysis are scale invariant from the Abell 
radius to about 20 hr 1 Mpc (Xu, Fang and Wu 1998, hereafter XFW). We will show that the scale 
dependence of M/L ratio is completely determined by the scaling factors of RF and MF of these 
multi-scale identified objects. If the mass field traced by the observed galaxies can match with 
popular models of CDM cosmogony, one can calculate the M/L's scale dependence by comparing 
the observed RF scaling with the model-predicted MF scaling. Using the data of the IRAS 1.2 Jy 
galaxy survey, we find that the M/L ratio is approximately flat on scales up to, at least, 24 h l 
Mpc. 

In §2 we calculate the mass functions of multi-scale identified clusters from simulation and 
calculate their scaling indices. In §3 we derive the richness functions and their scaling indices of 
multi-scale identified clusters from IRAS 1.2 Jy galaxy survey. In §4 we study the scale dependence 
of mass-to-light ratio by matching the scaling behaviors of theoretical mass functions with those 
of observational richness functions. 55 is our conclusion. 



2. Mass function of r^-clusters from simulation samples 

The model samples of mass distribution used in this paper are similar to XFW. They are 
generated by N-body simulations with the P 3 M code developed by Y.P. Jing (Jing & Fang 
1994). Three popular models of the CDM cosmogony, i.e. the standard CDM (SCDM), the 
flat low-density CDM (LCDM) models and the open CDM model (OCDM), are employed. The 
Hubble constant h (= i?o/(100 km/sec/Mpc)), mass density flo, cosmological constant Q,\ and 
the normalization of the power spectrum of the density perturbation ag are chosen to be (0.5, 
1.0, 0.0, 0.62) for SCDM, (0.75, 0.3, 0.7, 1.0) for LCDM, and (0.75, 0.3, 0.0, 1.0) for OCDM. The 
simulation parameters are: 1.) box size L = 310 h^ 1 Mpc; 2.) number of simulation particles 
N p = 64 3 ; and 3.) the effective force resolution rj = 0.24 h^ 1 Mpc. In this case, each particle 
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has a mass of 3.14 x 10 13 f2M^ 1 M & , which is small enough to resolve reliably objects with mass 
M > 5.5 x 10 14 /i _1 M Q .We produce 8 realizations for each model. 

For the simulation results, the clusters on scale r c i and with mass larger than M are identified 
by the multi-resolution (or multi-scale) analysis based on discrete wavelet transformation (DWT) 
(Fang & Thews 1998). It has been shown that the DWT-identified clusters on the Abell scale, i.e. 
r a fe-clusters, are the same as those given by friends- of -friends (FOF) identification (XFW; Jing & 
Fang 1994). That is, in terms of the r a b-cluster identification, the DWT method is as good as the 
FOF. Thus, the DWT method can be reliably used to systematically identify objects on scales 
r d > r ab (Xu, Fang & Deng 1999). 

With the multi-scale identification of simulation samples, we have model-predicted the mass 
function of r c /-clusters, n(> M,r c i), which is the number density of r c ;-clusters with masses larger 
than M. The MFs of r d -clusters with r d = 1.5, 3, 6, 12 and 24 h~ l Mpc in models of the SCDM, 
LCDM and OCDM are shown in Fig. 1. The middle part of mass function, roughly in the range 
1(T 4 > n(> M,r cl ) > 1(T 7 h 3 Mpc" 3 , is more reliable, as all the Poisson errors are relatively 
small. 

It can be seen from Fig. 1 that all the MFs have a similar shape. This is because of the MF 
scale invariance (XFW) as 

n(> M,r cl ) = n(> XM,r ab ), (1) 

where the scaling factor \{r c i/r ab ) is only a function of the ratio r c i/r ab , but independent of M. 
Obviously A = 1 when r c ; = r ab . The scaling factor A is an important statistic description of 
clustering on large scales. 

3. Richness function of r c ;-clusters from IRAS sample 

We apply the DWT multi-resolution identification to the redshift survey of the Infrared 
Astronomical Satellite (IRAS) galaxies with flux limit of f§o > 1.2Jy (Fisher et al. 1995). The 
IRAS sample is of 2 (physical space) + 1 (redshift space) dimension. A DWT algorithm for 
identifying r c ;-clusters from the 2+1 dimension samples has been developed (XFW; Xu, Fang 
&; Deng 1999). To minimize the radial effects, the data are divided into redshift intervals of 
[2500,5000], [5000,7500] and [7500,10000] km s _1 . We then identify r c/ -clusters in each interval 
on scales r c \ = 1.5, 3, 6, 12, and 24 h~ l Mpc. 

The sampling rate of the IRAS sample is low. The IRAS clusters identified on Abell scale r ab 
contain much fewer galaxies than optical clusters, and IRAS galaxies do not populate the clusters 
in the same manner as optical galaxies (Fisher et al. 1995; Webster, Lahav & Fisher 1997). This 
gives rise to higher Poisson errors. However, if we are interested in clusters larger than the Abell 
scale, this problem would be softened as structures contain more member galaxies. Meanwhile, 
IRAS galaxies are mostly of late-type. It traces clusters on various scales r c i uniformly. Therefore, 
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the samples of IRAS r c ;-clusters may be less affected by morphological segregation, and good for 
studying the scale-dependence of structures. 

With the identified r c ;-clusters of IRAS galaxies, one can find the richness function of IRAS 
r c ;-clusters, Nj R (> N g , r c i, z), which is the number density of r c /-clusters containing > N g members 
of IRAS galaxies in the redshift shell z. The RFs of r^-clusters with r c \ = 1.5, 3, 6, 12 and 24 h^ 1 
Mpc are shown in Fig. 2, in which the errors are Poissonian (Xu, Fang & Deng 1999). Despite the 
errors are large, it is clear that the shapes of all the RFs are basically the same. This feature can 
be described as a RF scaling invariance in each z shell, 

N m (> N g ,r d , z) = N IR {> X IR N g , r ab , z). (2) 

The scaling factor X[ R (r c i/r ab ) is a function of the scale ratio r c i/r ab , but independent of N g . 
Obviously Xj R = 1 when r c ; = r ab . The scaling factor Xj R of IRAS r c /-clusters as a function of 
fell fab is presented in Fig. 3. It shows that Xi R is also independent of z within the 1-a error bars. 



4. Scale dependence of mass-to-light ratio 

By definition, the mass-to-luminosity ratio for IRAS r c ;-clusters is the ratio between mass 
and luminosity of r c /-clusters. One can introduce the M/L ratio of the IRAS r c ;-clusters by a 
mass-number relation as 

M = R(r d )N G , (3) 

where M is the total mass of a r^-cluster, and N G the total number of member galaxies of the 
r c ;-cluster, defined by 

N G = N g /4>{z), (4) 

where (j>(z) is the selection function (Fisher et al. 1995). The definition of R{r c i) by eq.(3) has 
implicitly assumed that R(r c i), or the M/L ratio, is independent of the mass of the objects, since 
the M/L ratio of r a fe-clusters is found to be mass-independent on average (Bahcall & Cen 1993). 
Such independence may be caused by the fact that the scale of luminosity-related hydro dynamical 
process are smaller than that of clusters. Thus, luminosity functions and M/L ratios are 
statistically independent of environment on cluster scale. Therefore, it would be reasonable to 
assume that the mass-independence of the M/L ratio still hold on scales larger than clusters. 

The number-mass conversion coefficient R(r c i) is proportional to the mass-to-light ratio 
M/L(r c i) for IRAS galaxies on scale r c /. We have then 

M/L(r d ) = R(r cl ) 

M/L(r ab ) R(r ab Y [ > 



Using eqs. (3) and (4) , one can find the mass function of IRAS r c /-clusters njf>(> M,r d ,z) 

as 

n IR {> M,r ch z) = N IR (> {M<t>{z) / ' R(r cl )),r cll z). (6) 
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Therefore, the RF scaling (2) will further yield a MF scaling of IRAS r c /-clusters as 

n IR (> M, r d , z) = n IR (> Xf R F M, r ab , z), (7) 
where the MF scaling factor Xf R F is given by 

Xf R F (r cl /r ab ) = A M (r d /r ai )M, (8 ) 

or 

M/L(r cl ) = Xf F {r cl /r ab ) 
M/L(r ab ) Xi R {r d /r ab ) ' 

Namely, the scale-dependence of the M/L ratio is completely determined by the ratio between the 
MF and RF scaling factors of the IRAS r c /-clusters. 

However, despite that we have the observed scaling factor X[ R (r c i/r ab ), we are in lack 
of the observed Xf R \r c i/r ab ). This problem can be overcome by matching the IRAS galaxy 
traced mass field with popular models of CDM cosmogony. With properly selected parameters, 
CDM cosmogony models have been found to provide good approximation to many observational 
properties of large scale structures, especially the mass function of optical and X-ray clusters (e.g. 
Jing & Fang 1994; Bahcall, Fan & Cen 1997, and references therein). On the other hand, some 
individual clusters and superclusters from optical surveys have been identified from the IRAS data 
(Fisher et al. 1995; Webster et al. 1997). Therefore, it is reasonable to assume that the mass field 
of CDM cosmogony model provides a description for the mass field traced by the IRAS galaxies. 
Thus, the statistical properties of the simulated r c ;-clusters should be the same as that of the 
IRAS r a b-clusters. We have then 

X{r cl /r ab ) = Xf R F (r cl /r ab ), (10) 

i.e. the scaling factors of n(> M,r c i) and nj R (> M,r c i) are the same. Using eqs. (9) and (10) we 
have 

M/L(r cl ) = X(r cl /r ab ) 
M/L{r ab ) X IR {r cl /r ab y 

Thus, the M/L(r c i) ratio can be determined by the ratio between the model-predicted X{r c i/r ab ) 
and the observed Xi R (r c i/r ab ). 

The ratio X/X[ R as a function of r c i/r ab is plotted in Fig. 4. For LCDM and OCDM, the 
value of this ratio does not depart from 1 with a significance larger than 1-a in the entire range of 
scales from 3 to 24 h- 1 Mpc. In the case of SCDM, only at r d = 24 h 1 Mpc, the ratio shows a 
difference from 1 with 2 a level. Particularly, in the redshift shell [2500,5000] km s _1 , the ratio is 
very flat for all the three models. Thus, for IRAS sample, the M/L ratio is basically flat on large 
scales from 3 to 24 /i -1 Mpc, i.e. 

M/L(r cl ) 



M/L{r ab ) 



1. (12) 
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This result is weakly dependent on specific parameters of models of the CDM cosmogony. The x 2 
goodness is larger than 98% for all the three CDM models. 

Nevertheless, in Fig. 4, we can see a slight, but systematic increase of M/L with scale. A 
linear fitting by log(^r (r c i)/^(r a i,)) ~ a + 61ogr c / yields x 2 goodness larger than 99.9% for all the 
three models. The least-square fitted slope b is found to be 0.19 ± 0.17 for OCDM, 0.21 ± 0.17 for 
LCDM, and 0.41 ± 0.17 for SCDM. That is, the M/L ratio on scale ~ 20/1" 1 Mpc is probably 
larger than that of r a i, by a factor about 1.7 for LCDM and OCDM. (or, a factor of 2.9 for 
SCDM.) These results for OCDM and LCDM are consistent with the M/L ratio measurement 
of the Corona Borealis supercluster. Namely, M/L scale-dependence given by IRAS galaxies is 
consistent with optical data. 

5. Conclusion 

The RF scaling factor of multi-scale identified r c /-clusters is fundamentally important as a 
statistical measure for clustering on scales larger than clusters. Since these RFs depend only on 
the number of member galaxies of r c /-clusters, they can be directly detected from a galaxy redshift 
survey regardless the mass determination of individual objects. 

The M/L ratio of gravitationally clustering systems which consist of luminous objects and 
dark matter depends only on the scale of the system but not on total mass (Bahcall and Cen 1993). 
This is actually a scaling of cosmic gravitationally clustering. Therefore, the scale dependence of 
M/L ratio is closely related to other scaling indices of these objects. We have shown that the scale 
behavior of M/L is determined by the RF and MF scaling factors of r c ;-clusters. 

By comparing the mass function scaling of simulated samples with the richness function 
scaling of IRAS 1.2 Jy galaxies, we conclude that the assumption of M/L-flatness is basically 
correct in the range from 1.5 to 24 h^ 1 Mpc. This conclusion is weakly dependent on the 
cosmological parameters we have considered. 

We would like to thank an anonymous referee for a detailed report which improves the 
presentation of the paper. 
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Fig. 1. — Mass functions of clusters identified with radii r d = 1.5, 3.0, 6.0, 12 and 24 h Mpc for 
models of SCDM, LCDM and OCDM at the epoch z = 0. n(> M,r c {) is the number density of 
r c /-clusters with masses larger than M in unit of /i 3 Mpc -3 . M is in unit of h~ l M . 



Fig. 2. — Richness functions of multi-scale identified r c ;-clusters from IRAS galaxies in the redshift 
shells of [2500,5000], [5000,7500] and [7500,10000] km s _1 , with r d = 1.5,3.0,6.0 and Uh' 1 Mpc. 
N g is the number of member galaxies of the IRAS r c ;-clusters. nm{> N g , r d ) is the number density 
of r c /-clusters with number of member galaxies more than N g . 



Fig. 3. — Scaling factor \iR,(r c i/r a b) as a function of r c \. The error bars come from the clusters 
identified from IRAS galaxies. The triangle, circle and cross are, respectively for the redshift shells 
[2500,5000], [5000,7500] and [7500,10000] km s -1 . The triangle and cross data are horizontally 
shifted, respectively, by -O.Oldex and + O.Oldex for clarity. 



Fig. 4 — The ratio \{r d /r ab )/\ IR (r d /r ah ) or (M / L{r d )) / (M / L{r ab )) as a function of scale r d . 
The symbols are the same as those in Fig.3. 
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